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It is shown and discussed how the conventional system of hydrodynamic equations for ferrofiuids 
was derived. The set consists of the equation of fluid motion, the Maxwell equations, and the magne- 
tization equation. The latter was recently revised by Felderhof [Phys. Rev. E 62, 3848 (2000)]. His 
phenomenological magnetization equation looks rather like corresponding Shliomis' equation, but 
leads to wrong consequences for the dependence of ferrofluid viscosity and magnetization relaxation 
time on magnetic field. 
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I. INTRODUCTION. EQUATION OF FLUID 
MOTION 



o-ik = -pSik + V ( + ^ j + - I^ik), (1) 



In a recent paper [1], Felderhof made an attempt to 
revise the conventional hydrodynamic equations for fer- 
rofiuids. He proposed some modification in the equation 
of ferrofluid magnetization. A complete set of fcrrohy- 
drodynamic equations was worded for the first time in 
our paper [2] almost 30 years ago. In the Comment, we 
will consider the object of Felderhof's criticism, analyze 
his proposition, and explain why is it wrong. 

The main peculiarity of ferrofiuids is a specific relation 
between the magnetic and rotational degrees of freedom 
of suspended magnetic grains of which the fluids are com- 
posed. Therefore the concept of internal rotation first 
applied to ferrofiuids in [2] has proved to be very fruit- 
ful. The model [2] takes into account that the volume 
density of the angular momentum of ferrofiuids consists 
of both the visible ( "orbital" ) and the internal ( "spin" ) 
parts. The former, L = p(r x v), is associated with the 
translational motion of magnetic grains and molecules of 
the solvent. The latter, S, is caused by the rotation of the 
grains themselves and should be treated as an indepen- 
dent variable along with the fiuid velocity v, density p, 
and pressure p. However, an appropriate thermodynamic 
coordinate is the difference S — Ifl where ft — ^ curl v 
is the local angular velocity of the fiuid and / means the 
volume density of the particles moment of inertia. For a 
suspension of spherical particles / = ps4>(P /lO where 4> is 
the volume fraction of the dispersed phase, ps the parti- 
cles material density, and d the mean particle diameter. 
In this case it is convenient to set S — lujp where ijjp 
is the macroscopic (i.e., averaged over physically small 
volume) angular velocity of the particles. Any deviation 
of ijjp from gives rise to dissipation processes due to 
redistribution of angular momentum between L and S 
forms. (The angular momentum conservation law refers, 
sure, to the total angular momentum L -I- S ). These pro- 
cesses contribute the stress tensor cr^fc. For an ordinary 
(nonmagnetic) suspension the tensor has been derived by 
the methods of irreversible thermodynamics in [3] : 



where 

Sik = eikiSi , flik = (dvk/dxi ~ dv.Jdxk)/2 = e^^i , 

and tiki stands for antisymmetric unit tensor. Apart from 
the viscosity 77 , Eq. (1) contains once more kinetic coef- 
ficient: the spin relaxation time Tg = I/Grjip = psCp/GOrj . 
For d = 10 nm and rj = 10~^ Ps this formula gives 
Ts ~ 10"^^ s(!). Thus, the difference u:p — ft instantly 
decays whereupon the hydrodynamic description is re- 
duced to the common set of hydrodynamic equations. 
Ferrofiuids, however, give us an opportunity to maintain 
this difference by an extraneous magnetic torque which 
acts directly upon the particles rotation: 



6T](j){u}p - rj) = Mx H. 



(2) 



Here H is the magnetic field within the fiuid and M is the 
ferrofiuid magnetization. At the equilibrium in a station- 
ary field, M is described well by the Langevin formula 

Mo = nmL(0^, C = ^, L{C) = coth^ - r\ (3) 

where m is the magnetic moment of a single particle and 
n the number density of the particles. Eliminating the 
last term in (1) with the aid of the torque balance equa- 
tion (2) and including in aik the Maxwell tensor of mag- 
netic field, one gets [2,4] 

(T^k = -p5^k + ^(9^+5^) + 2^^'^^Hk - MkH,) 

+ ^{H,Bk-]-H^S,k). (4) 
47r 2 

On substitution Bk = Hk + AuMk in this tensor, we are 
convinced of its symmetry. Equation (4) and the momen- 
tum conservation law 



dvj _ duik !L - ^ ^( T7\ 

^ dt ~ dxk ' dt ^ dt^^ 
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determine the equation of ferrofluid motion 

= - Vp + r^V^v + (MV)H+i curl (M x H) . (5) 

In the calculation of the divergence of the stress tensor 
we have used the equations 



divv = 



curlH = 



divB = 



(6) 



i.e., the ferrofluid is considered to be incompressible and 
non-conducting. 

The system of equations (5)-(6) is still not complete 
since it does not determine the ferrofluid magnetization. 
The latter influences the fluid motion (see (5)) and de- 
pends itself on the motion as well. There arc two ba- 
sic ways to derive the missing magnetization equation. 
Both the ways have been proposed by the author with 
co-workers [2,5] and discussed in reviews [4,6,7]. 



II. PHENOMENOLOGICAL MAGNETIZATION 
EQUATION 

Originally the magnetization equation has been de- 
rived phenomenologically [2] as a modification of the De- 
bye relaxation equation [8] . To get the generalized equa- 
tion, one should introduce a local reference frame S', in 
which the suspended particles are quiescent on the av- 
erage, i.e., u^'p = 0. It is natural to assume that the 
magnetization relaxation is described in the system by 
the simplest Debye-like equation 



_ = --(M-Mo) 



(7) 



with Mo from (3). Other words, it assumes that any 
deviation (either in direction or magnitude) of M from 
its equilibrium value Mq decays according to the sim- 
ple exponential law (M — Mo) ^ cxp(— t/r) . Here 
T = ^'qV/ksT stands for the Brownian time of rotational 
particle diffusion {V = nd^/Q) since the particles are as- 
sumed to be rigid magnetic dipoles whose reorientation 
is possible only with rotation of the particles themselves. 
The frame of reference S' rotates with respect to the 
fixed ("laboratory") system S with the angular velocity 
Up . The rates of change of any vector A in systems E 
and S' are related by the kinematic expression 



dA 

'dt 



: X A - 



d'A 

~dr 



(8) 



Substituting here A = M , a;p from (2), and d'M./dt 
from (7), we obtain the equation sought: 

^ = nxM-- (M-Mo) - -?--Mx(Mx H) (9) 

dt T Grjcp 

(Shliomis, 1972). The last (relaxation) term in this equa- 
tion describes a process of approach of the vector M to 



its equilibrium orientation without change of the length 
of this vector. Equations (5), (6) and (9) constitute the 
complete set of conventional ferrohydrodynamic equa- 
tions. 

Let us compare Eq. (9) with the equation [1] 

^ = J2xM-7H(H,-H)-g^Mx(MxH) (10) 

(Felderhof, 2000). Here 'Jh is a positive phenomeno- 
logical constant and "H( = MC(M) is expressed in 
terms of the local magnetization" [1], that requires an 
introduction of a very inconvenient notation C{M) = 
{UbT /mM)L~^{M /nm), where L~^{x) means the func- 
tion inverse to Langevin function (3). On the face of 
it, Eq. (10) is only much less convenient than (9) since 
the left- and right-hand sides of (10) contain different 
relaxing values: M and H;(M). Felderhof claims how- 
ever that his "relaxation equation was derived from irre- 
versible thermodynamics (IT) , and differs from that pos- 
tulated by Shliomis. The two relaxation equations lead 
to a different dependence of viscosity on magnetic field." 
Let us consider both these statements. 

From the point of view of IT, the relaxation term 
(Mg — M)/t in (9) is neither more nor less "postulated" 
than the term 7//(H — H;) in (10). It is worth to remind, 
there are two methods of IT, and the both was proposed 
by Landau. The first of two, LI, based on conservation 
laws and the condition of positive entropy production, 
was applied for the first time to the building of hydrody- 
namics of helium [9], then hydrodynamics of fluids with 
internal rotation [3], liquid paramagnets [10], and some 
other liquids. Relaxation equation (10) has been also de- 
rived by the method. One has to be skilful enough to use 
this cumbrous method because it does not lead to the 
one and only form of sought equations. In this sense the 
second method, L2, is much more definite, simple, and 
direct than LI. It was first applied just to the description 
of relaxation of the order parameter in a non-equilibrium 
system [11]. An equilibrium value of the parameter (Mq 
in our case) corresponds to the minimum of an appro- 
priate thermodynamic potential $ (usually the Gibbs or 
Helmholtz free energy) depending on the magnetization 
M and other thermodynamic variables. Thus, at the 
equilibrium 9i>/9M = 0. Out of equilibrium this con- 
dition is not satisfied, so the relaxation process occurs: 
M changes in time approaching Mq. For small devia- 
tions from equilibrium, the derivative 9$/3M and the 
relaxation rate dM./dt are small. The relation between 
the two derivatives in the Landau theory is reduced to 
simple proportionality: 



dM 
"dt 



-1 



9$ 
dM. 



(11) 



with a constant coefficient 7 > 0. Hence we have 



d^ a$ dm /d^\^ „ 

— = = -7(^) <0 (12) 



dt dM dt 
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as it should be: when a system moves to equilibrium, 
its free energy decreases. In the case of a weakly non- 
equilibrium state of the system, one can substitute in 
(11) and (12) the expansion 



■Mo) 



where siibscript marks the point of equilibrium. As 
the first derivative in this point is equal to zero and the 
second one is positive, Eq. (11) turns into Eq. (7) with 
r-i = 7(52$/SM2)o and Eq. (12) takes the form 



dt 



(M - Mo)^ 



(13) 



Thus, Eq. (7) and hence Eq. (9) are well corroborated 
by the method of IT. Equation (10) does also not conflict 
with IT. Nevertheless it is wrong. As we show below, it 
leads to anomalous result for ferrofluid viscosity and mag- 
netization relaxation time. The pitfall of IT is discussed 
in Sec. IV. 

The Einstein formula for viscosity of suspension, = 
770(1 + 2.5(/)), was obtained without taking into account 
the rotational motion of suspended particles relative to 
carrier liquid. If however the particles angular velocity 
does not coincide with the angular velocity of the 
fluid n , there arise friction forces which manifest them- 
selves in an additional (so-called rotational) viscosity rjr ■ 
As the difference Wp — is maintained by the magnetic 
torque (see (2)), rotational viscosity turns out to be a 
function of the dimensionless field strength ^. In a sta- 
tionary field, the steady solution of Eq. (9) yields in the 
linear order in Or 



M - Mo = T±(nxMo) , Tj_ = 



2t 



2 + ei(0 ' 



(14) 



where r_L is the relaxation time of the transverse (to the 
field) component of the magnetization. For the Poiseuille 
fiow or the planar Couette flow under the field directed 
along the flow (i.e., H ± ), we find 



Mx H= -rj_Moi?ri, 



(15) 



while for arbitrary orientation of magnetic field the right- 
hand side of the expression should be multiplied by sin^ a 
where a is the angle between vectors H and . Let us 
substitute the magnetic torque (15) in (5). When the 
field H is homogeneous, the magnetic force (MV)H in 
the right-hand side of Eq. (5) vanishes, while two other 
terms of the equation may be grouped: 

'nW\+\ curl (M X H) = (77 -h i T±MoH) V^v . 

The quantity added here to the ordinary viscosity should 
be regarded as rotational viscosity 



Substituting here Mq from (3) and t± from (14), we de- 
rive the formula [2] 



^'-(^^=2'^'^2Tii(0 



3 ^-tanh^ 
= 2^'^fTTi^- ^''^ 

In the absence of magnetic field an individual particle 

" rolls" freely along corresponding shear surface with an- 
gular velocity u)p equal to , so that rir{0) = 0. Con- 
versely, r7r(0 attains its limiting value 



r]r{oo) = 3r](f)/2 



(18) 



(the saturation) when rolling of the particle is replaced by 
slipping: the field of sufficiently large intensity guarantees 
constancy of the particle's orientation, not allowing it to 
twist with the fluid. Note that the result (18) does not 
depend on a concrete form of the magnetization equation 
but follows directly from the equation of fluid motion (5). 
Actually, in the limit under consideration tjjp = , so that 
Eq. (2) takes the form M x H = -67?(/)J7 . Substituting 
this torque in (5), we immediately arrive at (18). In- 
deed, the value was obtained by Hall and Busenberg [12] 
as early as 1969 without use of any magnetization equa- 
tion. In any case, however, such an equation must not 
contradict the saturation value (18). Our formula (17) 
does satisfy the Hmit (18), whereas Felderhof's equation 
(10) leads to the quite difii'erent result. His final formula 
in [1] gives the value 



p 3 (nm)'^ 
Vr (00) = 7;V(f> 7—7 -7 ^ 



(19) 



which is evidently less than the correct value (18). Ac- 
cording to [1], 7jj = x/t where x = nmP/SksT is the 
initial magnetic susceptibility. Substituting the jh and 
(f) = nV in (19) we find GrjcfyjH = 2{nm)^ /3 , afterwards 
Eq. (19) yields r/f (oc) = 9t/0/1O . Thus the ratio of the 
Felderhof's limiting value of viscosity to the correct value 
(18) is equal to 0.6. Other words, in the limit ^ = 00 
Felderhof's equation (10) predicts Wp = 0.4 O (?!) in- 



stead of < 



as it must be. 



Vr = jT_lMoH. 



(16) 



III. MAGNETIZATION EQUATION DERIVED 
MICROSCOPICALLY 

Both the above-mentioned phenomenological methods 
allow to obtain linear relaxation terms in hydrodynamic 
equations [like (Mq — M)/r in (9)] and corresponding 
quadratic terms for the rate of the entropy growth or the 
free energy diminution (like that in (13)). It is clear that 
such terms are valid only for small departures from equi- 
librium. Indeed, Eq. (9) describes well the rotational 
viscosity for arbitrary intensity of a stationary magnetic 
field but small values of f2r (see, e.g., a good agreement 
between McTague's experiment [13] and Shliomis' the- 
ory [2]), or for small dimensionless amplitude ^ or fre- 
quency CUT of an alternating magnetic field (see, e.g., ex- 
periments on the negative viscosity and their explanation 
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in [14]). Meanwhile, to describe successfully the nega- 
tive ferrofluid viscosity at finite values of the parameters, 
we did need to use in [15,16] a more precise magnetiza- 
tion equation. Such a macroscopic equation should be 
derived from the kinetic Fokker-Planck equation which 
provides the microscopic description of particle diffusion 
in colloids. The program was realized by Martsenyuk, 
Raikher and Shliomis [5] soon after the phenomenologi- 
cal magnetization equation (9) was derived in [2]. 

The Fokker-Planck equation for a ferrofluid moving in 
a field H has the form [7] 

dW - " 
2r— =R(R-2Tf2-ex^)W, (20) 

where e = m/m is the unite vector along a particle mag- 
netic moment, ^ = mH/fcsT , and R = exd/de is the in- 
finitesimal rotation operator. Equation (20) determines 
the orientational distribution function W{e,t) of parti- 
cles magnetic moments. The macroscopic magnetization 
is determined by the relation M(t) = nm{e) where angu- 
lar brackets denote statistical averaging with the distri- 
bution function. Multiplying Eq. (20) by e and integrate 
over the angles, we arrive at the equation 



2t 



d{el 
dt 



2TOx(e)-2(e)-(ex(ex^)), (21) 



which however is not closed. Indeed, along with the first 
moment of the distribution function, (e), Eq. (21) con- 
tains the second moment (the last term in the equation). 
It is easy to make sure that the equation for the second 
moment includes the third one, and so on, thus there is 
the infinite chain of cross-linked equations. Ideally, how- 
ever, one would like to have only one equation since only 
the first moment magnetization - has a clear physi- 
cal meaning. An original scheme of closure of the first- 
moment equation (21), titled the effective field method, 
has been proposed in [5]. Let us explain the fruitful phys- 
ical idea. 

In equilibrium (fi = 0) under a constant magnetic; field 
the stationary solution of Eq. (20) is the Gibbs distribu- 
tion 



Wo{e) = 



47rsinh^ 



exp(^e) . 



(22) 



An averaging of the vector e with function (22) gives 
expression (3) for the equilibrium magnetization: 



Mo = nmL(0^/^ 



(23) 



Only in true equilibrium the magnetization is one or an- 
other function of the field. In a non-equilibrium state 
there is no connection between M and H : any arbitrary 
magnetization may be created - in principle - even in 
the absence of the field. Nevertheless, one may consider 
any value of M as an equilibrium magnetization in a cer- 
tain - specially prepared - magnetic field. This effective 



field He is related to the n on- equilibrium magnetization 
by the equilibrium relation: 



M = nmL(ee) 



(24) 



During the equilibrium settling process, the dimension- 
less effective field ~ mll,,/kBT tends to the true field 
^, so the magnetization (24) relaxes to its equilibrium 
value (23). Comparing (23) and (24), we see that the 
latter is obtained by averaging of e with the distribution 
function 



T^e(e) 



47r sinh 



exp(Cee) : 



(25) 



which differs from the Gibbs distribution (22) by replace- 
ment of the true field by the effective one. Carrying out 
the averaging in (21) with the function (25), we find the 
sought equation [5] 



d_ 
di 



'I 



Tie 



€eX(^eXO- (26) 



This equation defines the dependence of the effective field 
upon time, true field ^ , and the fluid vorticity il. Its 
solution being substituted into (24) determines the 
magnetization of a moving fluid. In the case of small 
departures from equilibrium, the effective field might be 
represented as a sum of the true field and some small 
correlation: = ^ + f. Then from (23) and (24) in the 
linear approximation in u we get 



M - Mo = nm[L'{£)u\\ 
where the components 



(27) 



are parallel and perpendicular to the true field, respec- 
tively. Employing the relation (27), one can reduce Eq. 
(26) to the linear magnetization equation 



dM 

dt 



f^xM- 



H[H(M-Mo)] Hx(MxH) 



(28) 



where relaxation times of the components of magnetiza- 
tion are 



dlnL(C) 

= — T 



Substituting t± from (29) in (16), we obtain [5] 



3 ^AlLiU 



(29) 



(30) 



In the same approximation our phenomenological equa- 
tion (9) also takes the form (28) but with other relaxation 
times: tm = r and t±_ is defined in (14). 
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FIG. 1. Reduced rotational viscosity 2r]r{^)/Srj<f> as a func- 
tion of ^ given by the effective field method [5] (Eq. (30), solid 
curve), by the phenomenological approach [2] (Eq. (17), dot- 
ted curve), and by Felderhof's approximation [1] (Eq. (32), 
dash-dotted curve). 

Figure 1 shows that though at first sight functions (17) 
and (30) do not appear alike, they agree fairly closely in 
the entire range of their argument. Both they are in 
a good agreement with experimental data of many au- 
thors and with computational results provided by direct 
numerical integration of the Fokker -Planck equation in 
linear approximation in fir [17]. Quite recently Felder- 
hof [18] solved this linearized equation by the Galerkin 
method using a large number of trial functions (the asso- 
ciated Legendre functions) . Comparing his "exact result" 
for the rotational viscosity with (30) and (17), he wrote 
that "the result of Martsenyuk, Raikher and Shliomis [5] 
is quite a good approximation, but the result of Shliomis 
[2] deviates up to 17 percent" . It is a pity that Felderhof 
forgot to mention in [18] his own prediction [1] - see the 
bottom curve in Fig. 1 - which deviates up to 40 percent! 



IV. DISCUSSION 

It is easy to see that Felderhof's "local field" H; in 
(10) is none other than our effective field He deter- 
mined by (24). Hence one can use the relationship (27) 
for a linear analysis of Eq. (10). Substituting in (27) 

= - ^ = {m/kBT){Ue - H), we find 



(He-H)ii 



(M-Mo)ii 



(He)j 



3xi(0 



(31) 



Let us substitute these relations in Felderhof's equation 
(10) and put there = x/t '■ with such a choice his 
equation coincides with (9) in the limit ^ -C 1 as it should 
be. Then for arbitrary ^ we obtain from (10), making use 

of (31) and (16), 



6L{0t 



e [2 + 3X2(0]' 



=2^-^2 + 3X^(0 ■ 



(32) 



Both the relaxation time and viscosity are wrong. In- 
deed, in a strong field, ^ ^ 1 , they take the magni- 
tudes = 6t/5^ and i]^ = 9rj(p/W , while it should be 
r_L = 2r/^ = Q'qV/mH and rjr = 2>ri(f)/2. The depen- 
dence J?^(0 shown in Fig. 1 strongly differs from two 
other curves. 

Let us give consideration to the question, why does 
Felderhof's equation lead to the anomalous results. We 
have shown above how does Debye equation (7) originate 
from the potential $(M). One can choose, however, as 
an independent variable the effective field and introduce 
the potential <I>(He). Then instead of (11) we obtain in 
similar fashion 



dt 



7 > 0. 



Acting further by the L2 method, we arrive at the equa- 
tion (cf. (7)) 



dt 



- (He - H) , 

r 



(33) 



where we set 7 ^ = (92<&/9He)T. Under the choice, Eq. 
(33) turns into (7) in the low field limit. 

Equations (7) and (33) satisfy the principal proposi- 
tions of the theory of linear response according to which 
the time rate x of change of a value x at each a mo- 
ment is determined by the value x at the same moment: 
X = x{x). Then, if x weakly deviates from its equilibrium 
value xq, one can expand x{x) over x and confine oneself 
to the linear term: x ~ —\{x — xq), where A is a positive 
constant. Thus, t in (7) and (33) should be considered 
as a constant. This inference seems important because 
the method LI of irreversible thermodynamics does not 
allow in principle to determine the field dependence of 
kinetic coefficients such as 7^ in (10). In the rotating 
reference frame E' the equation reads 



d'M 
dt 



-7H (He - H) . 



In contrast with (7) and (33), this equation relates one 
value (He) with time rate of change of another one (M). 
Therefore, under the nonlinear magnetization law (24), 
the coefficient 7^ wittingly cannot be constant but rep- 
resents a unknown function of ^ . 

Equation (33) has been written out in a coordinate 
system E'. Reverting to the immobile system E by the 
general formula (8) and eliminating u>p with the aid of 
(2), we obtain [19] 



rfHe 
dt 



= J2xHe--(He-H 



1 

&ri4> 



HeX(MxH). 



(34) 



This equation determines together with (24) the magneti- 
zation M in an implicit form, effective field He being the 
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parameter. In the case of small departures from equilib- 
rium, Eq. (34) can be linearized with respect to He — H 
and M — Mq. Using relationships (31) and 



dM 

IF 



dM dHe 



3X 



d(He)|, , L{0 (i(He)j 



dt 



+ 



dt 



we turn to Eq. (28) with t\\ = r and t±_ from (14). Thus, 
in linear approximation Eqs. (9) and (34) coincide with 
each other. As the result, both the equations yield the 
same relationship (17) for the rotational viscosity of fer- 
rofluids. 

It is worth to note, if Felderhof had used He as an 
independent variable correctly^ he would have arrived at 
Eq. (34). However, he has missed the opportunity. 



V. CONCLUSION 

Thus, Shliomis' theory consists of hydrodynamic and 
Maxwell equations (5)-(6) plus a magnetization equation. 
There are three kinds of the latter: (9), (26), and (34). 
It is well-established that Eq. (26) derived by the effec- 
tive field method from the Fokker-Planck equation yields 
quite accurate results for real ferrofluids. Indeed, a direct 
numerical simulation of the magnetic moment Brownian 
dynamics performed by Cebers [20,21] in the middle of 
80th has indicated that Eq. (26) describes perfectly the 
fluid magnetization in a wide range of parameters ^ and 
VLt. The same conclusion has been made in [6] under 
comparing of the solution of (26) with the results of nu- 
merical integration of the non-stationary Fokker-Planck 
equation (20). At the same time, the calculations [20,6] 
have shown that the phenomenological ccjuation (9) is 
valid for any field magnitudes ^ but only small enough 
fluid vorticities: Or < 1. Hence Eq. (9) can be recom- 
mended to the description of weakly noncquilibrium situ- 
ations, as the equation is far simpler for analysis than Eq. 
(26). The latter, however, guarantees the correct quanti- 
tative description of magnetization processes even if de- 
viations from the state of equilibrium are large, Or ^ 1, 
that is when Eq. (9) leads to erroneous results. Inter- 
estingly, our latest calculations have shown that the new 
equation for the effective field (34) is free from such a 
shortcoming: it is valid even far from equilibrium. There- 
fore, taking into account that (34) is nevertheless simpler 
than (26), one should recommend Eq. (34) to the most 
wide applications. 



As for the Felderhof 's equation (10), it does not stand 
up to comparison even with our phenomenological equa- 
tion (9) to say nothing of the microscopically derived Eq. 
(26). We have shown that incorrectly derived Eq. (10) 
leads to anomalous results (19) and (32) - see also Fig. 
1 - and that is why it should be rejected. 
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